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The invariant information introduced by Brukner and 
ZeiUnger, Phys. Rev. Lett. 83, 3354 (1999), is reconsidered 
from the point of view of quantum state estimation. We show 
that it is directly related to the mean error of the standard 
reconstruction from the measurement of a complete set of mu- 
tually complementary observables. We give its generalization 
in terms of the Fisher information. Provided that the opti- 
mum reconstruction is adopted, the corresponding quantity 
loses its invariant character. 

Invariants are important concepts of physics, because 
referring to them, physicists build up theoretical mod- 
els making it possible to explain consistently their ob- 
servations. Recently a project has been undertaken to 
build up the kinematics of the quantum theory from the 
information-theoretical principles • As a part of the 
project a new measure of the classical information gained 
from a measurement on a quantum system was intro- 
duced. This quantity summed over a complete set of 
mutually complementary observables (MCO) exhibits in- 
variance with respect to unitary transformations applied 
to the state of the system and/or to the measured set of 
MCO. Moreover, when properly normalized, it also quan- 
tifies the (maximum) information content, and therefore 
evaluates the processing power of physical systems. 

A nice feature of the invariant information of Brukner 
and Zeilinger is that its definition is operational. It is 
obtained by synthesizing the errors of a specially cho- 
sen set of measurements performed on the system. In 
this contribution, we analyze the invariant information 
from a different perspective. If certain observations are 
made on the system the obtained results can be in a 
natural way put together to form our estimate of the 
quantum state of the system. This hints on the ex- 
istence of a tight link between the information gained 
from a particular set of measurements and the error of 
the reconstruction based on the obtained results. Being 
motivated by the estimation theory we will show how 
to synthesize information gained from individual mea- 
surements in more general situations, namely, when (i) a 
complete set of MCO is not available, and when (ii) the 
observables measured are not necessarily mutually com- 
plementary (still non-commutative but not "maximally" 
non-commutative). Special attention will be paid to the 
invariance properties of the total information gained from 
a measurement. As a by-product, we will show an alter- 
native interpretation of the invariant information and will 
discuss the role that MCO play in the state estimation. 

Let us consider a measurement of an observable A = 
with a non-degenerated spectrum acting in 



the Hilbert space of dimension p. Each outcome j is 
detected with probability pj = TrpIIj. Brukner and 
Zeilinger's lack of information associated with such a 
measurement is defined as a sum of variances of indi- 
vidual outcomes per particle Ea = J2j'^j/-^^ where N 
is the number of particles incident on the measurement 
apparatus ||^. The total lack of information about the 
system is then obtained as a sum of those measures over 
a complete set of MCO, 
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where a = 1 . . . p -|- 1 and j = 1 . . .p label complemen- 
tary observables and their eigenvectors, respectively. The 
invariant information is nothing else than a properly nor- 
malized complement of the total lack of information E . 
Assume that a complete set of MCO exists for the given 
dimension p, which means there are p{p + 1) projectors 
Ua] satisfying (|J|l 
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One easily finds that the total error has already the de- 
sired invariance under the choice of the complete set of 
mutually complementary measurements, and under uni- 
tary transformation of the true state p of the system: 



E = p — Tip 
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The main argument for just summing up the errors of 
complementary observations is that MCO are indepen- 
dent in the sense that a measurement of one of them gives 
no information about the rest [Q. This is the standard 
interpretation of MCO observables based on the following 
statistical reasoning: Equation (^) could be interpreted 
as the conditional probability p{a, j\(3, k) corresponding 
to the detection of the quantum variable {a,j} provided 
that {/3, k} was true. This is symmetrical with respect 
to the interchange of observation and conditioning. Con- 
sequently, if a particular observable a is measured, the 
Bayesian posterior distributions associated with the ob- 
servables P a complementary to it are flat. Such ob- 
servation seems to enhance our knowledge just about a 
particular detected output without adding anything new 
to the unobserved variables. For this reason, MCO are 
sometimes called "unbiased observables" |8|. As will be 
shown in the following such a classical interpretation of 
MCO variables should be modified. 

There is a close link between the total error E and the 
error of the standard quantum tomography. Instead of 
discussing the errors of measured projections separately. 
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one can make a synthesis of all the results by forming 
an estimate of the unknown quantum state, and then 
evaluate the average error. Let us apply this strategy 
to the detection of a complete set of MCO. A generic 
quantum state may be decomposed in the basis of MCO 
as follows 



(4) 



where Ilaj are projectors on the eigenspaces of a complete 
set of MCO, the coefficients Waj being determined by the 
true probabilities paj = TrlpXIaj} as follows |^, 
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Now suppose that the complementary observations are 
done with N particles each. Due to fluctuations, the reg- 
istered relative frequencies faj will generally differ from 
the true probabilities Paj and so the estimated state 
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Given the same data, the accuracy of our estimate 
based on them strongly depends on the chosen recon- 
struction procedure. Direct inversion (|^) is simple and 
straightforward because it implicitly uses the apparent 
statistical independence of MCO. But this is not the only 
possibility there. Keeping in mind than we want to char- 
acterize the information gained through the measurement 
one should use the best reconstruction method available. 
To evaluate the error of the optimal estimation proce- 
dure it is more convenient to decompose the true density 
matrix in a basis of orthogonal observables rather than 
MCO, 



(11) 



where the — 1 unknown parameters Ok provide a min- 
imal representation of the state p, and Tj are orthonor- 
mal: TrlF^T/j} = Sj^- When the measurement is over an 
estimate p of the true state is formed. 



P = - + V'afc(/fc)rA 
P , 
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Let us see how much. Perhaps the most simple estima- 
tion strategy is the direct inversion based on the above 
mentioned independence of MCO, 
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The error will be quantified by evaluating the Hilbert- 
Schmidt distance between the true and the estimated 
state 



Tr{(p-p)n. 
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The mean distance (error) is then given by averaging d 
over many repetitions of the estimation procedure, each 
yielding slightly different estimates of Waj , 



Scst = (rf). 
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Using Eqs. 



(^) and the definition of complementarity 
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in Eq. (| )^ we get 
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Notice that for the given number of input particles the 
mean distance of the estimated state from the true state 
is proportional to the total error (|l|), which depends on 
p only through Trp^. This is one of the key results of 
this contribution. It shows an alternative interpretation 
of the total error E and the invariant information: The 
total lack of the information as defined by Brukner and 
Zeilinger determines the mean error of the standard re- 
construction based on the measurement of a complete set 
of MCO. 



The parameters specifying the estimated state depend 
on the registered frequencies according to the given re- 
construction procedure. Now we calculate the mean dis- 
tance between the estimated and true states. Trivial cal- 
culation shows that it is given by a sum of the variances 
of the estimated parameters. 



{d)=J2{{Aa,f 
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To proceed, the variances appearing in ( p3| ) must some- 
how be determined. This may be a nontrivial problem 
because generally the estimated parameters Uk might de- 
pend on the measured frequencies fk in a very compli- 
cated way. However, for our purpose, it is enough to eval- 
uate the performance of the optimum estimation. It was 
shown that the variances of estimated parameters cannot 
be less than the Cramer- Rao lower bound jl^. Further, 
it is known that maximum-likelihood estimators attain 
the bound asymptotically (for large N) |l^ . In our case 
the Cramer-Rao bound reads: 



((Aafc)') > F 
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where F is the Fisher information matrix defined as 

US 



Fki = 
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— logP(n|a)— logP(n|a) 
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Here P(n|a) denotes the probability of registering data 
n provided the true state is a, and the averaging is done 
over the registered data. For large N, one is allowed to 
replace the multinomial distribution P by its Gaussian 
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approximation. Keeping the notation of MCO registra- 
tions we have, 



logP(n|a) « -N^J2 
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where a^j = Npaj{l — Paj) as before. Under this ap- 
proximation the Fisher matrix becomes, 



Fki = N' 
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The main formal result of the paper is obtained from 
Eqs. (^3|) and (|l^): Operational information defined as 
the mean error of the optimal estimation from the mea- 
surement of the chosen set of observables is given by the 
trace of the inverse of the Fisher matrix, 
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It quantifies the maximum amount of information about 
an unknown state gained by the measurement performed. 
Eq. (^8|) has the following simple geometrical interpreta- 
tion: On registering counts n, the probability (^6|), con- 
sidered now as a function of a, characterizes the likeli- 
hood of various states a. In terms of the Fisher matrix 
it reads: 
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where a specifies the Maximum-likelihood solution. Let 
us define the error volume as the set of density matrices 
whose likelihoods do not drop below a certain threshold, 
C > const. According to (19) the error volume is an ellip- 
soid, the lengths of its axes being inversely proportional 
to the eigenvalues of the Fisher matrix. The optimal esti- 
mation error (d)opt can thus be interpreted as the sum of 
half-axes of the error ellipsoid which corresponds to the 
chosen set of measurement. From this point of view, the 
synthesis of all the quantum observations is equivalent to 
the registration of the orthogonal observables F' defining 
those axes. In this new representation the Fisher matrix 
attains the diagonal form, which means that the esti- 
mates of the transformed quantum-state "coordinates" 
a'f. = TrjpF'} fluctuate independently. Their statistical 
independence justifies the adding their variances as it ap- 
pears in Eq. (p^). 

This geometrical construction also shows that the 
structure of uncertainties related to the quantum tomog- 
raphy is richer than what a single number (|l^) might 
indicate. The chosen measurement might provide dif- 
ferent resolutions in different directions, depending on 
the shape and orientation of the noise ellipsoid. Let us 
mention in passing that the resolution optimized over all 
measurements has been shown to provide a natural sta- 
tistical distance between two quantum states |l^ . 

Now, we will go back to the measurement of a complete 
set of MCO. It is not difhcult to see that the invariance 



of the error (0) of the direct inversion (|^) stems from 
the fact that the projectors measured are considered as 
statistically independent observables. Each parameter ak 
in the decomposition (|^) is determined by measuring only 
one projector — the rest contributes nothing. As will be 
shown below, such data treatment cannot be optimal. 
This leads us to the following important question: Will 
the invariance of the estimation error survive if the direct 
inversion is replaced by the optimum data treatment? It 
turns out that the answer depends on the dimension p 
of the Hilbert space. In the simplest non trivial case of 
p = 2 the answer is in the affirmative. However for p > 3 
a deviation appears. Let us have a look at the behavior 
of (d)opt for p — 3. Straightforward calculation of Tr F^^ 
in this case yields: 
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Notice that the first expression on the right-hand side is 
exactly the invariant total error (|l^). The second ter m 
then quantifies the improvement of the optimum estima- 
tion upon the simple linear inversion (0). One can easily 
check that this term is not an invariant quantity. So it 
turns out that the error of the optimum estimation from 
the measurement of a complete set of MCO is not invari- 
ant with respect to unitary transformations. This means 
that the maximum amount of information about the true 
state that can be gained from such measurements might 
differ even for true states of the same degree of purity! 
One pays for the optimality of the reconstruction proce- 
dure by loosing its invariant character (universality). Or 
conversely, one can have an invariant reconstruction at 
the expense of loosing precision. 

Let us now address the broader aspects of optimal in- 
formation processing discussed above. It seems that the 
non-invariant character of ((i)opt originates from the "in- 
homogeneity" of the measured quorum of MCO. Two 
projectors drawn from the quorum can be either "comple- 
mentary", Tr{nin2} = 1/p, or orthogonal, Tr{nin2} = 
0, to each other. Consequently, the mesh established by 
MCO is not as regular as it might have seemed at the be- 
ginning. Provided that this is ignored, the observations 
may be regarded as "independent" on average, but such 
treatment is obviously not optimal. Some additional in- 
formation can be gained provided those observations are 
treated as noise dependent observations fl^ j. 

The invariance of the optimum estimation error can 
be restored provided a more "homogeneous" quorum of 
observables is measured. This can be illustrated on the 
thought measurement of the following observables. 
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where the normalization and the factor of 1 /2 ensure that 
< IIj < 1, so that Hj are actually POVM elements 
which can be associated probabilities 
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and (Zj specify the true state as before. Because we now 
have TrjlljTfc} — 5jk, the Fisher matrix ( p7| ) assumes 
the diagonal form. A simple calculation then yields 




Notice that this expression resembles the one for the in- 
variant total error (]^. The mean error of orthogonal ob- 
servations is obtained by substituting (t| = Npj{l — pj) 
and Eq. ^) in Eq. (||), 

N{d)or-tho = + - - Trp2. (24) 

2 p 

This is, obviously, an invariant quantity. 

Finally, let us mention that optimal quantum recon- 
struction schemes provide invariant estimation errors. 
Such an optimal measurement would consist in measur- 
ing the "unknown" quantum state along its eigenvectors, 
pllj = Xjllj. It would obviously give mutually exclusive 
results due to their orthogonality — Maximum likelihood 
estimation would always coincide with a deterministic 
data inversion. This gives the minimal achievable value 
of the average error: NEopt — X]j '^ji^ ~ ^j) — ^~ Trp^. 
It should be clear that in the case of a generic quantum 
state whose diagonalizing basis is unknown and should 
be estimated together with {Xj}, the optimal informa- 
tion can only be attained asymptoticalh^ in the limit of 
the large signal to noise ratio N —^ oo |l^] . 

To conclude, we have shown that the invariant infor- 
mation introduced by Brukner and Zeilinger is related 
to the problem of the estimation of a quantum state. It 
quantifies how the estimated state differs on average from 
the true states in the sense of the Hilbert-Schmidt norm. 
It depends on the quality of the measurement and on the 
data treatment adopted. Usually, the measurement of a 
complete set of mutually complementary observables is 
assumed. In this case the measured state may be easily 
found via a linear inversion. The average distance be- 
tween the estimated and true state is given directly by 
the invariant information. 

The amount of extracted information may be increased 
provided that data are used in optimal way achieving 
the ultimate limit given by (the inverse of) the Fisher 
information. When doing this the unitary invariance is 
sacrificed in favor of the accuracy of estimation. Conse- 
quently, some "directions" can be identified easier than 
another. This seems to be just an example of the in- 
terplay between the universality and effectiveness of a 
particular method. 

Finally, we touched the the concept of mutually com- 
plementary observables. One advantage of such obser- 
vations is that the data can be more easily processed 
to get information about signal via the direct inversion. 
This is possible, unitary invariant but not efficient — 



one can do something better. The inefficiency of such 
an approach can be recognized in the structure of the 
Brukner-Zeilinger information. In the full spectrum of 
mutually complementary observables some projections 
commute but some do not. However, their variances are 
always added regardless of the commutativity or non- 
commutativity! This is not optimal. In the light of this, 
the often discussed issue of the existence or non-existence 
of a complete set of mutually complementary observables 
in some dimensions seems to be just an academic ques- 
tion without essential consequences for the problem of 
the state determination. 
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